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ABSTRACT 


EXPERIMENTAL  INVESTIGATION  OF  GAS  BEARINGS 
WITH  ULTRA-THIN  FILMS 


ARON  SERENY 


The  experimental  investigation  described  here 
involves  the  highly  accurate  measurement  of  bearing 
clearances  on  the  order  of  10  microinches  in  self-acting, 
pivoted, narrow-slider  gas  bearings.  The  experimental 
measurements  are  based  on  light  interferometry  using  a 
variable-wavelength  pulsed  dye  laser  and  a CW  HeNe  laser 
as  monochromatic  sources.  The  light  interference  in  the 
gas  bearing  is  obtained  by  flying  the  slider  on  a very 
precise,  optically  flat  quartz  disk  through  which  the 
light  beam  is  transmitted. 

The  combined  effect  of  high  Knudsen  numbers  and 
surface  irregularities  on  the  flying  height  of  narrow  gas 
bearings  is  observed  by  varying  the  load  on  the  bearing 
and  the  ambient  molecular  mean  free  path.  The  experimen- 
tally measured  bearing  clearances  are  compared  quantita- 
tively with  rather  accurate  theoretical  predictions  obtained 
by  numerical  solution  of  Reynolds  differential  equation 

/ 


vii. 


for  compressible  fluids  with  slip  boundary  conditions. 

The  result  of  this  study  indicates  that  as  clearances 
in  narrow  gas  bearings  get  progressively  smaller,  while 
the  Knudsen  number  increases  to  values  beyong  0.1,  the 
theoretical  model  fails  to  predict  the  bearing  behavior. 

It  is  also  argued  that  this  failure  is  because  of  the 
weakness  of  the  continuum  model. 


viii. 


Introduction 


I . 


An  experimental  apparatus  to  study  the  lubrication 
phenomena  associated  with  very  low  clearances  was  design- 
ed and  fabricated.  It  is  now  working  as  a slider  bearing 
tester  at  the  Columbia  University  Lubrication  Research 
Laboratory.  Fig.  1 shows  the  experimental  setup  in  its 
present  state. 

The  experimenta]  portion  of  this  investigation 
involves  gas  bearings  operating  at  very  low  clearances 
such  as  those  often  encountered  in  gas  bearing  gyroscopes, 
flying  heads  for  rotating  magnetic  memories,  and  in  any 
gas  bearing  during  start  up  and  low  speed  operation.  The 
slider  bearing  used  for  this  research  is  a type  of  read- 
write  flying  head  used  in  many  conventional  computer  disk 
memory  units  (Winchester  head  ) . This  head  is  supported 
by  two  self  acting  air  bearings. 

Utilization  of  an  air  bearing  to  support  a magnetic 
recording  head  was  first  introduced  through  the  I.B.M. 
Ramac  disk  file  in  1956  [ 8]t  The  experimental  investi- 
gation by  Brunner,  Barker,  Houghton  and  Osterlund  [ 4], 
established  the  effects  of  various  parameters  on  the 
performance  of  air  lubricated  slider  bearings  with 
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relatively  large  minimum  film  thickness  (400  microinches) . 
This  experimental  work  compared  well  with  theoretical 
results  obtained  by  digital  computer  solution  of  Reynolds 
equation,  an  equation  derived  from  the  Navier-Stokes 
equations  [ 1] , which  applies  to  low  Reynolds  number  flow 
phenomena,  where  viscous  forces  are  predominant  over 
inertia  forces  and  where  clearance  variations  are  slow. 

At  that  time,  because  of  the  high  clearances  involved  (on 
the  order  of  100  molecular  mean  free  paths) , no  effects 
caused  by  rarefaction  were  considered. 

The  continuous  efforts  to  achieve  higher  recording 
bit  density  in  magnetic  recording  reduced  the  head  to 
disk  spacing  from  300  microinches,  reported  by  [ 4],  to 
the  nominal  gap  of  today's  flying  heads  (20  microinches) 

[ 2] . This  reduction  in  the  flying  height  along  with 
some  additional  development  in  head  design  increased  the 
recording  bit  density  from  1000  bit/inch-50  track/inch 
to  6000  bit/inch-400  track/inch.  The  relation  between 
recording  bit  density  and  head  clearance  is  given  in 
[11]  and  in  [ 8]  . 

The  lower  bearing  clearances  required  the  extension 
of  the  classical  Reynolds  equation  into  the  slip  flow 
regime,  Burgdorfer  [10] . The  suggestion  that  the  range 
of  applicability  of  the  Navier-Stokes  equation  can  be 
extended  somewhat  by  the  retention  of  slip  boundary 
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conditions  was  initiated  by  Schaaf  [23] . With  lower 
bearing  clearances,  the  gas  film  thickness  of  bearings 
in  use  today  is  on  the  order  of  several  mean  free  paths. 
This  modified  Reynolds  equation  is  used  for  the  design  of 
gas  bearings  where  the  dimensionless  ratio  of  the  ambient 
molecular  mean  free  path  over  the  bearing  trailing-edge- 
clearance,  >.  /h  , (Knudsen  number)  is  between  0.01  and 

u 1 

0.1.  The  effect  of  the  mean  free  path  being  on 
the  order  of  the  clearance  in  a spiral  grooved  thrust 
bearing  was  studied  by  Hsing  and  Malinoski  [19].  The 
results  reported  by  [19]  reveal  that  the  effect  of  slip 
boundary  conditions  could  contribute  a substantial 
reduction  in  performance. 

Only  very  recently  was  the  first  experimental  data 
provided  by  Tseng  [13]  who  reports  measurements  of  load 
versus  spacing  characteristics  of  gas  bearings  with 
clearances  of  twenty  microinches.  According  to  Tseng 
[13],  experimental  data  are  in  better  agreement  with 
numerical  solutions  when  the  slip  at  the  boundary  is 
included.  In  the  above  mentioned  report,  consideration 
is  not  given  to  crown  height,  nor  is  it  given  to  surface 
roughness  effects.  In  Tseng's  experimental  work,  mean 
free  path  effects  could  not  be  observed  for  Knudsen 
numbers  on  the  order  of  0.5  because  of  the  wide  slider 
bearings  that  were  used.  With  a wide  slider  bearing. 
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the  region  of  high  pressure  dominates  and  that,  in  turn, 
causes  the  load  carried  to  be  less  sensitive  to  slip 
boundary  conditions.  The  perturbation  solution  of 
Reynolds  equation  with  slip  boundary  condition,  as 
described  in  Appendix  D,  indicates  the  slip  effect  for  a 
wide  bearing  to  be  of  the  order  1/A,  (A  = bearing  number) . 
It  is  interesting  to  note  that  elementary  theoretical 
studies  performed  by  Christensen  [22]  and  by  Rhow  and 
Elrod  [15]  predict  an  increase  in  slider  load  with 
increasing  roughness,  while  Tseng's  clearance  measure- 
ments are  systematically  below  the  predicted  values. 

A review  of  the  techniques  for  the  measurement  of 
air  bearing  separations  is  given  by  Lin  [ 3] . These 
techniques  can  be  listed  as  capacitive,  inductive,  and 
optical.  With  larger  size  bearings,  such  as  the  3/4  in. 
dia.  slider  of  Ref.  [ 4],  the  capacitance  probe  technique 
was  reported  to  be  suitable.  For  slightly  smaller  slider 
bearings  (.42"  x .35"),  the  unshielded  capacitan-'e  probe 
technique  was  used  for  determining  slider  flying 
characteristics  and  for  measuring  flying  heights  in  the 
40  to  200  microinch  range  [ 6]  . 

The  implanted  capacitor  microprobes  used  by  Briggs 
and  Her)<art  [ 6]  were  almost  twice  as  large  in  diameter 
as  the  slcates  of  the  slider  bearings  used  in  today's 
flying  heads. 

The  miniature  size  of  today's  recording  heads  make 
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the  implantation  of  capacitance  probes  extremely  difficult 
if  not  impossible.  The  optical  light  interference  method 
is  applicable  to  such  small  slider  bearings.  This  method, 
as  applied  to  air  gap  measurements  between  a moving 
magnetic  tape  and  a static  recording  head,  is  described 
in  references  [16]  and  [17].  One  advantage  of  the 
interferometric  technique  is  that  it  results  in  a complete 
mapping  of  the  gap  distribution  between  the  magnetic  tape  or 
disk  and  the  recording  head.  A white  light  interferometric 
technique  applied  to  the  measurement  of  very  thin  films 
(10  microinches)  was  reported  by  Lin  and  Sullivan  [ 5] . 
Today  it  is  well  accepted  that  the  optical  interference 
technique  is  the  most  accurate  method  for  the  measurement 
of  ultrathin  gas  bearing  separations  [ 3]. 

In  this  experimental  research,  in  order  to  eliminate 
the  domination  of  the  high  pressure  region  and  to  allow 
for  significant  enough  slip  effects,  a relatively  narrow 
slider  bearing  is  used  (length/width  = 11).  No  experi- 
mental studies  of  narrow^  slider  bearings  operating  in  the 
slip  regime  have  been  reported. 

The  experiments  demonstrate  the  failure  of  the 
modified  Reynolds  equation  to  predict  accurately  the 
bearing  behavior.  The  results  indicate  that  as  the 
clearances  in  narrow  gas  bearings  get  progressively 
smaller,  while  simultaneously  the  Knudsen  number 
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increases  to  values  beyond  0.1,  the  theoretical  model 
fails  to  predict  the  bearing  behavior.  Also,  accurate 
measurement  of  the  crown  of  the  so  called  "flat"  slider 
and  the  effect  of  this  crown  height  is  given. 

In  addition  to  the  experimental  work,  attention  was 
given  to  improve,  and  get  some  further  insight  into,  the 
analytical  treatments  of  the  governing  lubrication 
equation.  Various  numerical  methods  were  investigated 
and  comparisons  were  made  with  exact  solutions, 
developed  here,  for  wide  slider  bearings.  This  has  been 
the  first  time  that  the  exact  solution  to  Reynolds 
equation  with  the  molecular  mean  free  path  effects 
included  were  carried  through.  Also,  in  order  to 
study  analytically  the  slip  phenomena,  an  asymptotic 
solution  was  obtained  and  it  indicates  effects  to  be  of 
the  order  1/A  for  the  wide  bearing  case. 
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II.  Theoretical  Background  of  Self-Acting  Slider  Bearing 
with  Ultra-Thin  Film 

The  theoretical  treatment  of  self-acting  gas  lubri- 
cated slider  bearings  is  well  established.  Thorough 
treatment  of  the  fundamentals  of  gas  lubricated  films  is 
given  by  Fuller  [ 1]  and  by  Gross  [18] . The  steady  state 
Reynolds  equation  is  given  here  in  dimensionless  form. 
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1 3r,  9P] 
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3Y 

a: 

Cl 

^1 

where  P = P/Pg*  H = h/hj^^^,  X = x/l,  Y = y/l, 

A = 6\iUl/  (Pgh^niin^  ' number 

= ambient  pressure;  h . = minimum  clearance 

mxn 

Its  validity  in  predicting  the  behavior  of  hydro- 
dynamic  bearings  with  continuum  laminar  flow  between  two 
smooth  surfaces  has  been  proven  repeatedly.  Equation  ( 1) 
describes  the  fluid  dynamics  of  bearings  with  clearances 
much  larger  than  the  molecular  mean  free  path  and  the 
size  of  surface  asperities.  Further  analysis  of  gas 
bearings  with  thinner  and  thinner  fluid  films  requires 
review  of  the  concept  of  continuum  in  the  description  of 
fluid  flow.  The  continuum  hypothesis  for  the  flow 
between  two  bearing  surfaces  weakens  with  clearances 
approaching  the  value  of  the  molecular  mean  free  path. 
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Gas  dynamics  is  usually  divided  into  the  regions  of 
continuum,  transition  and  molecular  flow.  In  order  to 
extend  the  range  of  validity  of  the  Reynolds  equation  to 
the  transition  region,  Burgdofer  [10]  introduced  the 
commonly  used  first  approximation  of  slip  at  the  wall. 
This  consists  of  using  the  new  boundary  condition 


U 


wall 


f (A,) 

a 


^ J wa  1 1 


(2) 


where  f{Aa)  = A-, 

U = velocity  tangent  to  the  wall 
y = coordinate  normal  to  the  wall 


o = reflection  coefficient,  [29] 


With  0=1  for  diffuse  reflection,  this  led  to  the 


following  modified  Reynolds  equation 


9 

9 X 


H’p(1  + 

' PH'9  X 


9 PH 
9 X 


(3) 


with  the  Knudsen  number,  k being  defined  as  the  ratio  of 
the  ambient  mean  free  path  over  the  bearing  trailing- 
edge-clearance . 

This  modified  Reynolds  equation  predicts  a reduction 
in  the  slider  load  carrying  capacity.  This  reduction  is 
pronounced  at  low  bearing  numbers.  Although  this  equation 
has  been  used  for  several  years  in  actual  design  and 
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analysis  of  low-flying  slider  bearings,  the  first 
experimental  study  by  Tseng  [13]  in  the  area  of  rare- 
faction came  only  recently  (1975).  Tseng  suggests  the 
introduction  of  a surface  accommodation  coefficient  of 
less  than  unity  in  order  to  explain  the  lower  load 
measured  by  his  experiment  in  comparison  to  the  numeri- 
cally predicted  one.  The  numerical  solution  of  the 
modified  Reynolds  equation  using  various  values  of 
molecular  mean  free  path  gives  the  load  versus  clearance 
behavior  described  in  Fig.  2.  This  figure  shows  a 
significant  reduction  in  minimum  clearance  under  constant 
load,  due  to  slip  at  the  boundary.  In  Fig.  2,  the 
Knudsen  number  increases  with  decreasing  values  of  trail- 
ing edge  clearance;  the  Knudsen  number  increases  from 
values  less  than  0.1  to  values  greater  than  0.3. 
Therefore,  since  the  Knudsen  number  is  larger  than  0.1, 
the  designer  should  expect  the  validity  of  these  design 
curves  to  weaken  at  clearances  below  20  microinches. 

In  addition  to  the  mean  free  path  effects,  the  very 
small  air  gaps  in  recent  usage  stimulated  the  initiation 
of  some  basic  studies  of  the  effects  on  bearing  perfor- 
mance due  to  the  roughness  of  the  stationary  and  moving 
parts.  Tseng  and  Saibel  [20],  Christensen  and  Tender 
[21] , and  Elrod  [14] , have  shown  using  Reynolds  equation, 
that  surface  roughness  can  considerably  affect  the 
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bearing  performance.  More  recently,  Rhow  and  Elrod  [15] 
showed  that  there  is  a general  increase  in  load  carrying 
capacity  due  to  striated  roughness  when  calculations  are 
based  on  mean  film  thickness.  They  carefully  distinguish 
between  calculations  based  on  mean  film  thickness  and 
ridge  film  thickness,  and  their  conclusion  indicates  that 
Stokes  typo  grooving  reduces  the  load  capacity  if  it  is 
calculated  on  the  basis  of  ridge  film  thickness.  If  the 
roughness  does  not  cause  excessive  rates  of  change  of 
film  thickness,  they  found  the  increase  in  load  carrying 
capacity  to  be  much  larger  for  bearings  with  the  rough- 
ness on  the  stationary  surface  than  for  the  case  of 
roughness  of  the  moving  surface  (Fig.  3) . No  studies 
have  yet  attempted  to  apply  theoretical  predictions  on 
roughness  effects  using  some  measured  input  from  an 
actual  bearing  surface. 

Consequently,  for  thin  film  bearings,  under  fixed  load, 
there  is  a decrease  in  clerance  due  to  approaching  rarefacted 
conditions  and  an  increase  in  clearance  due  to  surface 
irregularities.  Fig.  2 shows  the  predicted  flying  height 
without  including  surface  roughness.  Since  all  actual 
bearings  have  some  surface  irregularities,  the  curves 
in  Fig.  2 should  be  moved  toward  the  higher  load  - higher 
clearance  region.  Although  some  theoretical  treatments 
have  already  been  performed  on  certain  roughness  effects 
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as  well  as  on  slip  boundary  conditions,  as  yet  no  theore- 
tical studies  have  combined  rarefaction  and  roughness 
effects  or  predicted  the  regions  of  domination  for  gas  lub- 
ricated bearings.  In  this  experimental  study  the  observed 
behavior  of  a particular  slider  bearing  with  its  surface 
roughness  characteristics  is  given. 
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III.  Experimental  Apparatus 

In  order  to  test  the  performance  of  ultrathin  film 
slider  gas  bearings  with  various  mean  free  paths  and 
flying  heights,  a testing  setup  with  the  ability  of 
controlling  these  parameters  is  needed.  Since  the 
pressure  is  inversely  proportional  to  the  molecular  mean 
free  path,  the  ambient  pressure  in  the  testing  chamber  is 
to  be  held  at  sub-atmospheric  values  in  order  to  achieve 
the  desired  increase  in  the  mean  free  path.  To  control 
the  flying  height,  a variable  loading  mechanism  was 
constructed . 

The  testing  setup  was  designed  as  follows: 

* to  have  a vibration-free  flat  and  very  smooth  running 
surface  oposite  the  stationary  slider; 

* the  ambient  pressure  to  be  controlled  from  14.7  to 
1.5  psia ; 

* the  mean  free  path  to  be  controlled  from  2.6  to  26 
microinches ; 

* to  set  the  load  on  the  slider  bearing  from  10  to  100 


[ 


A 


grams ; 
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* to  have  the  very  clean  environment  needed  for  bearings 
with  ultrathin  films. 

A.  Description  of  Apparatus  and  Instruments 


The  experimental  setup  that  was  designed  and  built 
at  the  Columbia  University  Lubrication  Research  Labora- 
tory is  a highly  versatile  testing  bench  for  small-size, 
self-acting  gas  lubricated  slider  bearings.  The  essential 
components  of  the  apparatus  are  illustrated  schematically 
in  Fig.  4.  The  equipment  consists  of  two  major  parts: 
one,  the  controlled  environment  chamber  which  encloses 
the  running  and  stationary  parts  of  the  bearing;  and  the 
second,  the  optical  bench  which  carries  the  interferome- 
ter optics  and  the  image  recording  equipment  for  the 
measurement  of  bearing  clearances. 

The  chamber  itself  is  divided  into  two  enclosures  by 
a middle  plate.  This  middle  plate  acts  as  a partition 
between  the  enclosure  with  the  speed-controlled  DC  motor 
drive  and  the  enclosure  with  the  bearing  being  tested. 

By  enclosing  the  electric  motor  and  the  belt  drive  in  the 
environmental  chamber,  the  difficult  problem  of  sealing 
the  high  precision  spindle  is  solved. 

The  bearing  testing  apparatus  can  be  operated  with 
or  without  the  upper  box  depending  on  its  functional 


Pulsed  Monochromatic  Light  Source  i— — 

A)  Pulsed  Tunable  Dye  Laser  Recording  Film 
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Figure  4 Schematic  Diagram  of  Experimental  Apparatus 
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requirements.  For  operation  in  atmospheric  conditions, 
the  experiment  is  performed  without  the  upper  cover,  with 
the  testing  setup  positioned  in  front  of  the  laminar  flow 
air  bench.  For  operation  under  controlled  ambient 
conditions,  the  upper  box  is  positioned  so  that  it  covers 
the  bearing.  In  this  mode  of  operation,  the  necessary 
clean  environment  is  provided,  by  letting  the  gas  enter 
the  upper  chamber  through  a line  of  absolute  filters  and 
exit  from  the  lower  chamber  through  a vacuum  pump.  This 
cycle  provides  a continuous  flow  of  gas  from  the  upper  to 
the  lower  chamber  and  thus  prevents  contaminants  from  the 
belt  drive  and  the  electric  motor  from  entering  the  bear- 
ing test  area. 

In  the  sealed  environmental  chamber,  the  mean  free 
path  of  the  lubricant  can  be  varied  by  simply  changing 
the  pressure  and  the  lubricant  gas.  Because  flight 
behavior  at  high  Knudsen  numbers  is  of  interest,  the 
chamber  is  designed  for  subambient  operation  only.  In 
the  upper  chamber,  a highly  polished,  double  surface, 
optically  flat  disk  with  surface  flatness  of  X/20  is 
mounted  on  the  spindle.  The  optical  flat  is  made  of 
quartz  and  has  the  dimensions  of  1.75  inches  thick, 

10  inches  O.D.,  1.8  inches  I.D.,  and  it  weighs  10  lbs. 

The  disk  has  a peak  to  valley  surface  roughness  of  0.5 
microinches.  In  order  to  achieve  the  desired  vibration- 
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1 

free  operation  at  velocities  up  to  6000  R.P.M.,  the 
spindle,  with  the  optical  flat  at  one  end  and  the  pulley 
at  the  other,  was  dynamically  balanced  to  within  less 
than  100  microinches  peak  to  peak.  Fig.  5 illustrates 
the  balancing.  The  actual  balancing  was  performed 
using  24  balancing  screws  in  3 different  planes  along 
the  shaft  of  the  spindle.  The  electric  motor  with  the 
driving  pulley  was  also  balanced  to  within  100  microinches 
peak  to  peak  vibration  at  the  nominal  operating  speed. 

While  the  moving  part  of  the  bearing  is  the  quartz 
disk,  the  stationary  part  is  a spring  loaded  two-rail, 
tapered-flat  ferrite  slider  with  nominal  length  of  0.22 
inch  and  a width  of  0.02  inch  per  skate.  The  arm 
holding  the  gimbal  with  the  head  is  clamped  to  a precision 
stage  with  six  degrees  of  freedom:  three  rotations  and 

three  linear  displacements  around  the  x,  y,  and  z 
directions.  The  precision  stage  permits  the  alignment  of 
roll  and  pitch  angles  during  flying. 

The  variable  loading  mechanism  applied  to  the  slider 
consists  of  a pneumatic  rubber  diaphragm  with  a loading 
button.  The  load  on  the  slider  can  be  continuously  ad- 
justed by  pressure  variation  under  the  diaphragm.  The 
arm  is  made  of  Delrin*  and  is  instrumented  by  means  of 


Trade  name  of  polyformaldehyde,  Du  Pont  Co. 
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' Figure  5 View  of  Disk  Balancing 

four  strain  gages,  two  on  each  side,  which  measure  the 
load  with  a sensitivity  of  10.5  microinches/inch  strain 
per  gram. 

Hence,  an  addition  of  5 grams  to  the  load  causes  the 
I strain  reading  to  increase  by  53  microinches/inch.  The  arm 

and  the  loading  diaphragm  have  been  calibrated  by  means  of 
dead  weights.  The  loading  diaphragm  and  the  calibrated  arm 
are  illustrated  in  Fig.  7. 

In  this  experimental  setup,  the  bearing  film  thick- 
ness is  measured  using  optical  interferometry.  The 
usual  application  of  interferometry  requires  a transparent 
and  optically  flat  reference  plane,  a monochromatic  light 
source,  and  a reflecting  microfinish  on  the  surface  to  be 


A 


12  feet 


Figure  6 Schematic  Diagram  of  Interferometric  Setup 


Ferrite  Beryllium  Copper 

Slider  Gimbal 
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analyzed.  By  projecting  the  light  from  the  monochromatic 
source  through  the  flat  reference  plane  and  onto  the 
reflective  surface  to  be  analyzed,  interference  "fringes" 
between  the  incident  and  refracted  light  may  be  observed. 
The  position  and  number  of  fringes  is  a function  of  the 
distance  and  parallelism  between  the  two  surfaces.  Know- 
ing the  wavelength  of  the  monochromatic  light,  the  observ- 
ed surface  can  be  mathematically  described  in  numerical 
terms  as  a function  of  the  number  and  type  of  observed 
fringes . 

The  actual  interferometric  setup  is  illustrated  in 
Fig.  6.  The  .3  mwatt  CW  HeNe  laser  is  located  right 
behind  and  coll  inear  with  the  pulsed  tunable  laser.  By 
passing  the  continuous  monochromatic  source  through  the 
coaxial  flash  lamp  of  the  dye  laser,  the  CW  collimated 
light  beam  is  allowed  to  travel  through  the  same  path 
and  to  use  the  same  recording  setup  as  the  pulsed  laser. 
While  the  CW  allows  a continuous  observation  of  the 
interference  pattern,  th.e  pulsed  laser  records  only  the 
fringe  pattern  with  its  variation  in  a period  of  0.5 
microsecond.  The  blue,  4500  A and  the  red,  6328  A , 
wavelengths  give  the  distinct  difference  in  fringe  pattern 
needed  for  the  clearance  measurements.  Two  such 
interference  patterns  are  shown  in  Fig.  8. 
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B.  Description  of  the  Slider  Bearing  and  its 

Prcparat ion 
*-  - . 

The  slider  bearing  used  for  this  experimental  study 
is  a recording  flying  head  of  a conventional  computer 
disk  memory  unit.  A schematic  description  of  the  slider 
bearing  and  a drawing  to  actual  size  is  given  in  Fig.  9. 
The  design  concept  that  permitted  the  adaptation  of  the 
two-rail  tapered-flat  ferrite  slider  to  the  low-load  head 
was  developed  at  I.B.M.  [2].  The  most  significant 
design  difference  between  the  two-rail  slider  and  its 
predecessor  heads  is  in  the  suspended  mass  and  load  force. 
The  head  is  made  of  high  density  sintered  ferrite  with 
excellent  wear  and  surface  finish  characteristics, 
properties  highly  needed  to  accommodate  for  the  bearing 
to  start  and  stop  in  contact  and  to  maintain  a finite 
constant  spacing  between  the  head  and  the  recording 
surface.  The  head  dimensions  were  designed  such  that 
under  10  gram  load  the  actual  center  of  slider  pitch  and 
roll  is  at  the  recording  gap  that  is  located  between  the 
two  skates  at  the  trailing  edge. 

This  narrow  slider  bearing  with  a relatively  high 
length-to-width  ratio  allows  a significant  side  leakage 
of  lubricating  air  from  beneath  the  slider.  The  slight 
pressure  generated  by  the  tapered  portion  of  the  bearing 
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View  of  Typical  Interference  Patterns 
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Actual  Size 
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Figure  9 Dual  Skate  Tapered  Flat  Slider  Bearing 
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is  leaked  to  lower  values  before  being  compressed  again 
by  the  trailing  end  of  the  flat  portion  of  the  skates. 

The  integrated  pressure  profile  under  the  skates  results 
in  a load  of  10  grams  with  nominal  trailing  edge  clearance 
of  20  microinches  at  speeds  of  1500  in/sec. 

To  obtain  a qualitative  knowledge  of  the  bearing 
surface  roughness  characteristics,  a few  of  the  bearing 
surfaces  were  studied,  using  Scanning  Electron  Microscopes. 
Some  of  the  typical  micrographs  are  given  here  in  Figures 
10-11. 

Generally,  the  micrographs  exhibit  a surface  with 
two  types  of  irregularities:  large  ones  at  great 
intervals,  and  very  fine  ones  closely  spaced.  The  large 
grooves,  possibly  caused  by  crystal  pullouts  during 
lapping,  are  on  the  order  of  5-20  microinches  in  depth 
and  width,  spaced  at  intervals  of  200  microinches  or 
larger.  The  irregularities  spaced  at  small  intervals 
were  estimated  to  be  on  the  order  of  one  microinch  in 
depth  and  width. 

The  electrical  conductivity  of  the  ferrite  bearings 
is  sufficient  for  the  SEM  application  so  that  no  conduc- 
tive coating  was  necessary. 

One  of  the  major  challenges  in  testing  gas  bearings 


is  the  demanding  standards  of  cleanliness.  A study  of 
the  levels  of  surface  cleanliness  needed  when  assembling 


Figure  10  SEM  Micrograph  of  Slider  Surface  Figure  11  SEM  Micrograph  of  Slide 


29 


and  testing  gas  lubricated  bearings  is  given  by  Galvin, 
Monecraft  and  Patterson  [7] . The  methods  applied  in 
cleaning  the  skates  of  the  flying  head  for  the  present 
experimental  research  are  those  described  in  Ref.  [ 7]. 
The  cleaning  processes  include  organic  solvent  and 
deionized  water  washes  of  the  surfaces.  It  was  found 
necessary  to  repeat,  in  a trial  and  error  fashion,  this 
washing  process  while  a flying  attempt  was  performed 
between  each  wash.  During  a flying  try,  surface  contami- 
nation could  be  observed  on  the  projected  pattern  of 
interference.  The  appearance  and  accumulation  of  dirt 
between  the  initial  cleaning  processes  could  easily  be 
observed  during  the  flight.  It  should  also  be  reported 
that  if  no  dirt  marks  appeared  between  the  bearing 
surfaces  during  the  first  few  minutes  after  take-off, 
the  heads  flew  with  no  dirt-caused  disturbance  for 
several  hours. 
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C.  Accuracy  of  Measurements  and  Calibration  of 
Instruments 


A description  of  the  analytical  formulation  relevant 
to  this  experimental  study  has  been  given  in  Sections  II 

and  VI  of  this  report.  The  knowledge  of  the  following 

* 

quantities  is  required  for  the  correlation  of  the 
experimental  and  theoretical  results: 


(a)  u 
^a 

Pa 

^red 

''blue 


viscosity  of  air 

molecular  free  path  of  air  at  atmospheric 
pressure 

atmospheric  pressure 

wave  length  of  HeNe  monochromatic  light 
source 

wave  length  of  the  Pulsed  Dye  Laser  (coumarin) 


(b)  U 


h 

r 


linear  velocity  of  disk  under  the  slider 
slider  total  length;  (Fig.  9) 
tapered  length;  (Fig.  9) 

distance  of  pivot  point  from  trailing  edge 

width  of  skates;  (Fig.  9) 

crown  height  of  the  flat  portion  of  the 

bearing  sk.jtes;  (Fig.  9) 

ramp  height;  (Fig.  9) 
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The  values  of  the  quantities  in  list  (a)  were 

_ 8 

assumed  to  have  constant  values;*  m = 0.26  x 10  lb><sec/ 

, _ 6 . ^ 
in  , = 2.65  X 10  inches,  = 14.7  psia,  = 

6328  A,  and  X,  , = 4500  A.  All  these  values,  with  the 

blue 

exception  of  X.  , , are  believed  to  be  accurate  to  within 

^ blue 

1%.  The  wavelength  of  the  blue  monochromatic  light  source 
is  reported  by  the  supplier  (Phase-R,  New  Durham,  New 
Hampshire)  to  be  within  4%. 

The  values  of  the  physical  quantities  in  list  (b)  of 
this  section  were  determined  experimentally.  The  linear 
velocity  of  the  disk  under  the  slider  was  monitored  by  a 
magnetic  pick-up.  The  frequency  of  pulses  generated  in 
the  magnetic  coil  were  counted  using  an  electronic 
counter.  The  estimated  error  in  speed  measurement  was 
less  than  1%.  The  transistorized  closed-loop  control 
system  for  the  motor  speed  kept  the  set  speed  to  within 
a range  better  than  1%.  The  geometrical  dimensions  of 
the  slider  bearing  f,  s^,  s^  were  measured  with  an 

accuracy  of  .0005  inch,  which  gives  a maximum  relative 
error  of  3%.  These  estimates  are  based  on  applying, 
whenever  possible,  at  least  two  methods  of  measurement. 

The  ramp  height  was  measured  wiht  an  accuracy  of  ^j-g^/2  = 
12.4  (uin.),  which  gives  a maximum  relative  error  of  4%. 


* Ambient  temperatures  were  kept  in  the  range  of  70-74"F 
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The  crown  height  determination  suffered  from  the  largest 
relative  error  because  of  its  extremely  small  value; 

(3  microinches  1/10  of  1 wavelength  of  visible  light) 
the  measurements  of  this  crown  were  not  better  than  ±1.0 
microinches,  which  in  relative  error  terms  indicates  a 
possible  error  of  30%. 

The  load  on  the  slider  bearing  was  determined  by  the 
read-out  of  the  change  in  strain  in  the  Delrin  arm 
supporting  the  slider  bearing.  The  arm  was  calibrated 
while  loading  it  with  dead  weights.  To  minimize  drift  in 
the  strain  read-out  caused  by  the  slightest  temperature 
variation  related  to  the  low  thermal  conductivity  of 
Delrin,  gages  with  the  most  suitable  self-temperature- 
compensating  backing  material  were  chosen.  With  this 
proper  choice  of  backing  material,  the  drift  in  the  read- 
out and  the  error  in  the  calibration  curve  were  reduced 
to  a maximum  relative  error  of  less  than  10%  for  the 
lowest  load  of  10  grams.  The  instrumented  arm  was 
tested  for  variation  of  strain  because  of  air  drag  and 
forced  convection,  and  these  effects  were  included  in 
the  above  error  estimate. 

The  accuracy  of  the  flying  height  measurements  were 
of  the  order  *1.5  microinches.  For  the  larger  flying 
heights  this  error  is  only  of  order  10%  while  at  the 
lowest  trailing  edge  clearance  measurements  the  relative 
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error  was  about  30%. 
identification  three 


To  minimize 
methods  were 


errors  in 
carefully 


fringe-order 
followed . 


These  methods  are  described  in  Section  IV. 
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W . Experimental  Results 
A . General  Remarks 

I The  slider  bearing  used  for  this  experimental  study 

I is  a recording  flying  head  of  a conventional  computer 

* 

; disk  memory  unit.  The  heads  for  this  experimentation 

wore  provided  by  the  Advanced  Development  Laboratories 

of  Xerox  Corporation  and  C.D.C.  A schematic  description 

of  the  recording  head  and  a drawing  to  actual  size  are 

given  in  Fig.  9.  The  head  is  mounted  onto  a beryllium- 
I 

copper  gimbal  that  permits  the  head  to  incline  relatively 
freely  in  any  direction.  The  gimbal  is  welded  onto  a 
leaf  spring  that  loads  the  head  with  10  grams  through  a 
so-called  pivot  point.  This  10  gram  load,  provided  by 
the  spring,  is  the  load  under  which  all  take-offs  and 
landings  were  performed. 

After  the  static  measurements  of  the  slider  geometry 
(described  in  Sec.  IV-B) , the  gimbal-spr ing  load  was 
measured.  Attaching  the  gimbal-spring  to  the  loading 
arm,  the  head  was  lowered  against  a weight  measuring 
micro-balance  and  from  the  reading  on  the  instrument, 
the  stiffness  and  the  load  of  the  spring  were  established. 
Next,  the  loading  coll  and  the  instrumented  supporting 
arm  were  calibrated  against  dead  weigts.  It  was  found 
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that,  although  the  zero  reading  (no-load)  of  the  strain 
in  the  arm  varied  somewhat,  the  slope  of  the  load  versus 
strain  remained  within  tight  limits.  With  this  in  mind, 
the  zero  for  determination  of  the  total  applied  load  was 
determined  by  using  the  reading  under  the  ten  gram  load 
of  the  gimbal-spring . 

The  variation  in  strain  read-out  between  static  and 
dynamic  conditions  both  under  atmospheric  and  sub- 
atmospheric  bearing  ambient  pressure  was  investigated, 
and  it  revealed  changes  in  strain  of  ±5  microinch/inch 
(equivalent  to  i,5  gram).  The  strain  gage  instrumented 
arm  was  checked  for  performance  in  partial  vaccum  (up  to 
152  mmHg) ; no  variation  in  strain  read-out  was  observed 
due  to  the  change  in  ambient  pressure. 

Several  experimental  steps  were  performed  to  study 
in  general  the  bearing  behavior  in  this  setup.  The 
interference  pattern  under  static  condition  (disk  not 
spinning)  between  the  slider  and  the  disk  surface  was 
studied  for  take-off  clearance  configuration.  After  a 
few  initial  starts  and  stops  the  head  was  allowed  to 
land  on  the  slowly  coasting-down  disk.  The  white,  red 
and  blue  light  interference  patterns  were  applied  to 
investigate  the  relative  configuration  between  the  slider 
and  the  disk.  These  measuremnets  disclosed  no  visible  contac 
fringes  between  the  lOgr.  loaded  slider  and  the  quartz  disk 
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after  landing . The  clearance  at  the  leading  edge  of  the 
slider  was  'v.  10  microinches  while  at  the  trailing  edge 
on  the  order  of  2 microinches.  This  observation  of  no- 
contact take-off  was  confirmed  by  the  fact  that 
approximately  1000  take-offs  and  landings  performed 
during  this  experimental  work  and  no  visible  contact 
marks  on  the  two  sides  of  the  bearing  (slider  and  disk) 
were  produced. 

To  observe  stability  of  the  bearing,  the  quartz  disk 
and  the  air  gap,  a fast  pulse  (.45  microsecond)  dye  laser 
was  used.  The  laser  was  triggered  using  a signal 
generated  once  every  cycle  by  a magnetic  probe.  By 
delaying  the  triggering  signal,  pulsed  interference 
patterns  were  taken  at  different  points  along  the 
circumference  of  the  disk.  These  recorded  interference 
patterns  revealed  no  change  whatsoever  in  the  bearing 
clearance  topography.  This  indicated  an  insignificant 
level  of  vibration.  The  slider  followed  the  disk  with 
its  waviness,  if  any,  (flat  to  A/20)  extremely  well. 

In  the  load-vs-clearance  measurements  described  in 
Figs.  17  and  19  , for  = 14.7  psia,  two  sets  of 

di 

measurements  were  included.  In  each  of  these  figures 
one  set  of  the  measurements  was  done  without  the  upper 
cover,  while  the  other  set  was  done  with  the  cover  in 
place.  These  measurements  generally  show  no  variation 
in  performance  botwwen  the  covered  and  uncovered  testing. 
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B.  Measurement  of  Slider  Dimensions 

The  initial  step  taken  for  both  the  numerical 
solution  and  the  experimental  measurements  of  flying 
height  is  the  measuring  of  the  geometrical  dimensions 
of  the  slider  tested.  The  crown  height  in  the 
direction  along  the  skates,  the  pivot  location  and  the 
ramp  length  and  height  were  found  to  be  especially 
important.  The  dimensions  of  the  head  used  for  the 
experimentation  and  given  in  Table  1 were  measured 
using  optical  microscopes  and  interferometric  techniques 
with  a special  lOX  magnification.  The  head  itself  is 
described  in  Fig.  9. 

Except  for  the  ramp  height,  all  given  dimensions 
are  the  result  of  measurements  made  by  these  two  methods. 
The  dimensions  obtained  with  the  two  methods  were  in 
agreement  with  each  other  within  less  than  4%.  The 

i 

I 

dimensions  given  in  Table  1 for  the  slider  bearing 
tested  are  averages  and  round-offs  of  the  measurements. 

The  ramp  height,  one  of  the  most  important 
dimensions  for  bearing  performance  prediction  was 
measured  by  optical  means.  Leaning  the  slider  against 
an  optical  flat  while  reflecting  a monochromatic  light  at, 
and  back  from,  the  slider  surface,  brought  about  the 
highly  dense  interference  pattern  indicated  in  Fig.  12. 


i. 


Figure  12  Interference  Pattern  for  Ramu  Height 
and  Length  Measurements 
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This  interference  pattern  needed  at  least  lOX  magnifica- 
tion for  fringe  clarity. 

The  crown  of  the  skates  was  measured  using  a project- 
ed interference  pattern  with  lOX  magnification  obtained 
while  flying  the  head  at  a speed  of  1490  in. /sec.  loaded 
with  10  grams.  Fig.  8 shows  such  an  interference 
pattern.  The  gradual  increase  in  spacing  between 
consecutive  fringes  indicates  presence  of  surface  crown. 

At  this  point,  the  numerical  solution  for  the  pressure 
distribution  for  the  appropriate  speed  and  load  was  used 
to  check  for  possible  elastic  deformation  of  the  slider 
due  to  this  load.  The  various  calculations  with 
specific  assumption  of  beam  shape  and  load  distribution 
predicted  maximum  deflections  on  the  order  of  1/5 
microinch.  Additional  calculations  under  even  higher 
load  indicated  these  deformations  to  be  negligible. 

In  order  to  increase  the  number  of  interference 
fringes  within  the  skates  and,  by  that,  increase  the 
number  of  data  points  along  the  slider  and  the  accuracy 
of  crown  height  measurements,  a monochromatic  source  with 
the  shortest  wavelength  possible  is  needed.  The  optical 
alignment  and  image  recording  procedures  limit  the 
monochromatic  sources  to  be  in  the  visible  range  of  light. 
In  this  interferometric  setup  a short-pulse  tunable  dye 
laser  was  used.  A schematic  description  of  crown  height 


Figure  13  Crown  Height  Measurement,  U = 1490.  in/sec 
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measurement  is  given  in  Fig.  13. 

Some  crude  attempts  were  made  to  estimate  crown 
height  across  the  bearing.  Such  attempts  showed  this 
crown,  if  any,  to  be  on  the  order  of  1 microinch.  Numeri- 
cal solutions  including  this  1 microinch  crown  across  the 
bearing  did  not  show  significant  differences  in  perfor- 
mance. From  here  on  the  skates  were  taken  to  be  flat  in 
the  direction  perpendicular  to  motion. 

Finally,  the  lapped  slider  profile  used  for  all  the 
numerical  calculations  is  given  in  Fig.  14. 
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C.  Comparison  with  Numerical  Results  for  Various 
Bear inq  and  Kn udsen  Numbers 

This  section  describes  the  numerical  solution,  and 
the  experimental  technitjue;  a comparison  between  these 
results  is  given. 

The  numerical  solution  used  to  obtain  the  theoretical 
results  is  described  in  Appendix  A.  In  order  to  compare 
the  experimental  results  v/ith  the  theoretical  predicitons, 
an  accurate  solution  of  Reynolds  equation  with  slip 
boundary  conditions  was  needed.  Two  major  difficulties 
were  encountered  in  solving  Reynolds  equation  numerically 
for  this  particular  low-load  low-clearance  slider 
bearing.  These  problems  resulted  from  numerical  instabi- 
lities in  solving  Eq . ( 3)  with  large  bearing  numbers 

(A  'V  1000)  and  from  inaccuracies  in  load  predictions 
caused  by  the  discontinuities  in  clearance  slope. 

To  overcome  the  instabilities  in  the  numerical 
solution,  a careful  distribution  of  the  variable  grid 
spacing  was  needed.  The  l^nowledge  of  the  trailing  edge 
boundary  layer  thickness  being  on  the  order  of  1/A  was 
the  criterion  applied  for  the  positioning  of  the  grid 
points.  For  half  the  symmetric  bearing,  a mesh  grid 
15  31  was  used  (15  spacings  across  and  31  spacings 

along  the  bearing) . 
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To  handle  inaccuracies  in  pressure  profile  solutions 
caused  by  the  large  discontinuities  in  clearance  slope, 
the  integral  descretizat ion  described  in  Appendix  A was 
used.  This  algorithm,  obtained  by  using  the  Gauss 
integral  theorem,  satisfies  mass  balance  at  each  grid 
point  as  the  column  solution  method  is  applied  [ 9]. 

This  discretization  reduces  the  order  of  the  partial 
differential  equation  by  one.  With  this  discretization, 
positioning  of  grid  points  at  clearance  discontinuities 
eliminates  the  need  to  approximate  derivatives  at  lines 
of  singularity. 

The  numerical  solution  was  checked  against  the 
exact  1-D  solution  of  Reynolds  equation  with  slip 
boundary  condition  (Appendix  B) . The  numerical  program 
was  also  checked  against  results  obtained  by  a high-A 
algorithm  obtained  by  F.lrod  and  Cheng  [30]. 

Converrjonce  checks  using  larger  number  of  grid 
[loints  were  also  performed. 

These  tests  indicated  accuracies  of  2%  and  of  7% 
for  bearing  numbers  of  1000  and  4000,  respectively. 
Solutions  are  not  known  to  be  available  in  the  literature 
for  the  modified  Reynolds  equation  ( 3)  for  a narrow 
bearing  configuration. 

The  computer  program  using  IBM  FORTRAN-G  compiler 
needs  a core  of  160K  bytes  and  40  seconds  of  computation 
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time  for  each  solution  of  the  Reynolds  equation  for  the 
pivoted  slider  tested. 

The  objective  of  this  experimental  research  was  to 
compare  measured  bearing  clearances  with  theoretically 
predicted  values.  A flow  chart  describing  the  procedure 
for  the  testing  of  a slider  bearing  is  given  in  Fig.  15. 

After  the  cleaning  of  the  bearing  surfaces  and  the  taking 
of  the  static  measurements  (steps  which  are  described  in 
Secs.  III-B  and  IV-B)  the  head  was  flown  at  a high 
velocity  and  low  load  for  the  measurements  of  both 
slider  crown  and  slider  dimensions.  The  crown  measure- 
ment is  described  in  Sec.  IV-B. 

Next,  after  all  the  geometrical  dimensions  were 
determined,  a parametric  study  of  the  slider  performance 
was  conducted.  The  pivoted  slider  bearing,  whose 
dimensions  are  given  in  Table  1 and  in  Fig.  9,  was 

I 

; tested  under  various  ambient  pressures  and  molecular  mean 

J free  paths  with  loads  of  10  to  30  grams  and  disk 

velocities  of  800  to  2000  in. /sec. 

I One  of  the  most  important  steps  in  the  flow  chart 

I 

[ given  in  Fig.  15  is  the  processing  of  a picture  of  the 

* 

. interference  pattern.  Knowing  the  order  of  the  fringes 

appearing  in  a recorded  pattern  allows  for  clearance 
measurements  with  an  accuracy  of  2 microinches.  If  an 
error  is  made  in  the  fringe  order  determination,  the 

_ d 
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inaccuracies  in  clearance  measurements  could  go  to  values 
of  15  microinches.  Therefore,  determining  the  fringe 
order  has  major  importance  in  the  measurements  of 
clearance  topography. 

Three  independent  methods  were  applied  to  eliminate 
any  possible  error  in  fringe  order  determination. 

One  of  the  best  ways  to  determine  the  fringe  order 
was  to  follow  and  count  fringes  during  take-off,  loading, 
pressure  change  and  disk  speed  variation.  Being  able  to 
follow  the  fringes  during  these  transient  steps  was  one 
of  the  benefits  of  the  extremely  stable  and  flat  quartz 
disk.  A typical  stable  fringe  pattern  is  shown  in  Fig.  8. 
The  top  fringe  pattern  in  this  figure  includes  the 
effects  of  clearance  variation  with  time,  since  the 
exposure  time  is  on  the  order  of  a full  disk  rotation. 

In  previous  experimental  studies  by  Lin  and  Sullivan  [ 5] 
and  by  Tseng  [13]  on  the  performance  of  these  bearings 
at  very  low  clearances,  due  to  the  surface  waviness  of  a 
large  thin  disk  (14.0"  dia.,  1/4"  thick),  no  such 
stability  was  achieved.  The  disk  used  by  [ 5]  and  [13] 
represents  more  a functional  surface  than  an  optically 
flat  one. 

The  second  independent  way  to  determine  fringe 
order  took  advantage  of  the  fact  that  two  blue  fringes 
(4th  and  5th  blue)  will  appear  between  the  third  and  the 
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fourth  red  fringes.  The  order  in  which  the  blue  and  the 
red  fringes  line  up  is  given  in  the  left  column  of  Fig. 

16. 

The  third  method  applied  to  determining  fringe  order 
was  by  trial  and  error.*  Assuming  first  the  orders  of 
of  red  fringes  the  skate  profile  was  plotted  and  then  a 
check  of  this  profile  was  performed  using  the  blue 
fringes.  This  trial  and  error  was  repeated  two  or  three 
times  until  a fit  similar  to  the  one  in  Fig.  16  was 
achieved . 

Experimental  results  versus  numerical  predictions 
are  given  here  in  Figs.  17-25.  These  figures  describe 
the  bearing  load  (in  grams)  versus  the  trailing  edge 
clearance  (in  microinches)  for  various  ambient  pressures 
and  bearing  velocities.  On  the  upper  line  of  these 
plots  the  bearing  numbers  based  on  trailing  edge 
clearance  and  the  Knudsen  numbers,  k = X /h  are  given. 

The  ambient  molecular  mean  free  path  of  the  bearing 
was  set  by  lowering  the  ambient  pressure.  Assuming 
isothermal  conditions,  which  can  be  shown  to  be  valid, 
the  molecular  mean  free  path  varies  inversely 
with  the  pressure;  the  mean  free  path  was  set  to  values 

* Thi.<3  method  was  not  used  after  a while  since  the  first 
two  were  found  superior. 

*. 
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of  2.6,  5.2,  and  7.8  microinches.  With  trailing  edges 
on  the  order  of  5 to  12  microinches,  the  Knudsen  number, 
k , was  increased  to  values  beyond  one.  This  indicates 
that  load  versus  clearance  measurements  were  perfomed 
for  bearings  where  the  value  of  the  trailing  edge 
clearance  was  the  same  as  the  molecular  mean  free  path 
of  the  ambient  gas.  One  could  already  suspect  that  in 
low-load  bearings  the  theoretical  models,  based  on  slip 
flow  approximations,  will  not  perform  well,  especially  at 
the  trailing  edge  region. 

After  careful  observation  and  cross  plotting  of  the 
theoretically  predicted  and  experimentally  measured 
results,  several  observations  related  to  Figs.  17-25  can 
be  made.  First,  examine  Figs.  17  and  18,  for  results 
under  a constant  load  of  10  grams.  The  results  on  these 
figures  indicate  significantly  better  agreement  for  the 
bearing  with  the  2000  in. /sec.  than  for  the  800  in. /sec. 
velocity.  Note  that  the  bearing  with  the  higher  velocity 
runs  with  a clearance  almost  twice  as  large  as  the  one 
with  the  lower  speed.  In  Fig.  18,  the  measured  flying 
height  under  the  10  gram  load  is  in  relatively  good 
agreement  with  the  predicted  clearance.  This  measure- 
ment, which  agrees  the  best  with  theoretical  predictions 
already  has  a relatively  high  Knudsen  number  (k  = .11). 
The  experimental  setup  did  not  allow  for  bearing  testing 
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with  Knudsen  numbers  below  ,11, 

At  this  stage,  without  even  increasing  the  molecular 
mean  free  path,  results  indicate  that  as  the  clearances 
in  this  narrow  gas  bearing  get  smaller,  while  simulta- 
neously the  Knudsen  number  increases  from  ,1  to  ,173,  the 
theory  progressively  fails  to  predict  bearing  behavior. 
The  progressive  failure  of  the  theory  to  provide  the 
needed  side  leakage  to  lower  the  bearing  load  is  seen  in 
Fig,  18,  where  discrepancies  go  from  20%  to  25%  and  to 
40%,  while  the  Knudsen  numbers  vary  from  ,1  to  ,173  and 
to  ,210,  respectively.  Looking  now  at  Fig,  17,  where 
clearances  change  from  about  12  to  5 microinches,  this 
trend  of  progressively  increasing  disagreement  stops  at 
about  7 microinches  and  the  general  direction  seems  to 
reverse  itself.  Here,  with  clearances  on  the  order  of 
5 microinches,  the  agreement  seems  to  improve  as  the  gap 
gets  smaller  and  smaller.  At  this  point  of  the  analyis, 
thoughts  related  to  increasing  load  because  of  surface 
roughness  could  come  to  mind. 

In  Fig,  18,  an  attempt  was  made  to  use  the  suggested 
surface  accommodation  coefficient  by  [13]  but,  as  the 
figure  indicates,  it  did  not  change  the  results 
significantly  enough. 


In  Figs,  20-24  measurements  with  larger  molecular 
mean  free  path  are  given.  In  Figs,  30-32,  the 


i. 


53 


experimental  results  for  pressures  14.7,  7.35,  and  4.9 
psia  indicate  that  the  discrepancy  between  experiment  and 
theory  increases  when  the  ambient  pressure  is  reduced 
(i.e.  , ambient  mean  free  path  is  increased)  . In  Fig.  24  , 
the  bearing  reaches  bearing  and  Knudsen  numbers  on  the 
order  of  27000  and  1.5,  respectively.  This  discrepancy, 
under  low  load  at  a high  Knudsen  number,  could  also  be 
observed  in  Fig.  5 of  [13],  where  a wide  slider  bearing 
with  a very  low  load  was  tested. 

It  is  interesting  to  observe  that  slio  effects  are 
still  important  at  such  large  bearing  numbers  (A  2000)  . 
An  asymptotic  study  on  the  slip  effects  at  large  bearing 
numbers  is  given  in  Appendix  D.  This  analytical  solution 
indicates,  as  previous  researchers  show,  that  slip  effects 
diminish  at  large  values  of  A.  This  statement  immediately 
raises  the  question  of  how  large  is  a large  A.  For  a 
finite  width  bearing,  the  magnitude  of  A,  an  indicator 
for  slip  effects,  should  only  be  used  after  the  pressure 
profile  (or  load)  is  compared  with  the  solution  for 
Since,  for  the  narrow  slider  bearing  used  here,  the  load 
carried  is  only  a fraction  of  that  of  a A -►<»  bearing, 
bearing  numbers  on  the  order  of  2000  should  still  be 
considered  small. 

Finally,  a cross  plot  of  the  experimental  and 
theoretical  results  is  given  in  Fig.  25.  Curves  of  Fig. 
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25  represents  the  solution  of  Reynolds  equation  for 
various  values  of  bearing  and  Knudsen  numbers  for  a 
pivoted  slider  bearing.  This  figu -e  indicates  that,  for 
better  agreement  with  tne  experi  , the  slide  leakage 
in  the  bearing  ought  to  be  increased. 

Figures  26  and  27  present  the  numerical  solution  and 
the  experimentally  measured  pivot  clearance  versus  the 
load  on  the  slider.  These  two  figures  indicate  a good 
agreement  between  theoretical  predictions  and  the  experi- 
ment for  the  pivot  clearance. 
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Figure  23  Load  vs  Minimum  Clearance,  U = 1000  in/sec 
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Figure  26  Load  vs  Pivot  Clearance,  U = 1490  in/sec 


27  Load  vs  Pivot  Clearance,  U = 1490 
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V.  Remarks  and  Recommendations 


The  design  of  precise  gas  bearings  to  support  today's 
recording  heads  with  core  clearances  of  20  microinches 
has  been  greatly  aided  by  rather  accurate  theoretical 
predictions.  To  design  the  next  generation  of  heads 
flying  with  clearances  down  to  5 microinches,  the 
theoretical  models  must  be  modified.  The  following 
recommendations  are  given  to  outline  some  of  the 
additional  highly  needed  research  for  the  modelling  and 
design  of  ultra-thin  gas  bearings: 

* developing  ways  to  check  the  accuracy  of  side  leakage 
obtained  by  the  currently  available  numerical  methods; 

* developing  theoretical  modification  for  the  governing 
lubrication  equation  to  be  applicable  in  the  transition 
regime  of  gas  dynamics; 

* improving  and  furLher  developing  the  numerical  methods 
in  solving  Reynolds  equation  for  both  narrow  and  wide 
bearings  for  cases  with  bearing  numbers  reaching 
10,000; 

* performing  a theoretical  study  on  gas  bearings  with 
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both  slip  and  surface  roughness  effects  included  to 
observe  regions  of  domination; 

* bearing  surface  roughness  characterization  which 
provides  surface  data  adequate  for  the  theoretical 
modeling  of  roughness  effects  on  bearing  performance. 

In  conclusion,  the  author  believes  that  this 
investigation  has  contributed  toward  improved 
experimental  methods  and  has  provided  the  first  set  of 
experimental  data  on  low- load  low-flight  slider  bearings. 
It  is  hoped  that  this  experimental  analysis  may  serve  as 
the  first  link  in  future  theoretical,  numerical,  and 
experimental  work  in  the  design  and  development  of 
ultra-thin  gas  bearings. 
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APPENDIX  A 

Numerical  Solution  of  Reynolds  Equation  with 

SJ^ijD  Boundary  Condition  for  Cases  of  Large 
Bearing  Number  (\  > 300) 

The  numerical  solution  of  the  governing  steady  state 
lubrication  equation  that  has  been  applied  to  obtain  the 
theoretical  results  is  described  in  this  Appendix.  In 
order  to  compare  the  experimental  results  with  the  theo- 
retical predictions,  an  accurate  solution  of  Reynolds 
equation  with  slip  boundary  conditions  was  needed. 

Numerical  solution  methods  that  are  described  in  a 
review  paper  by  Castelli  and  Pirvics  [ 9)  were  applied  to 
the  modified  Reynolds  equation  for  slider  bearings  with 
large  bearing  number.  The  numerical  methods  described 
here  will  point  out  ways  to  overcome  major  difficulties 
caused  by  large  bearing  number,  A > 300,  and  because  of 
large  discontinuities  in  clearance  slope  of  a slider 
bearing.  It  should  be  noted  that  in  this  investigation, 
the  effort  to  obtain  an  accurate  numerical  solution  was 
carried  beyond  the  accuracy  needed  from  the  designer 
point  of  view. 

Generally,  there  are  two  distinct  methods  of 
discretization  applied  to  Reynolds  equation  ( 3) 
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Figure  lA  Differential  vs  Integral  Discretization 
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differential  discretization  and  integral  discretization. 
Each  method  applied  to  a differential  equation  of  the 
form  resembling  Reynolds  equation  is  described  in  Fig.  lA. 
The  integral  discretization  is  obtained  after  integrating 
the  Reynolds  equation  over  the  shaded  area  in  Fig.  2A  and 
then  applying  the  Gauss  theorem.  The  differential 
discretization  is  the  method  to  be  used  wlien  the 
differential  equation  is  written  in  a finite  difference 
form  in  a straightforward  fashion. 


Fig.  2 A 


The  governing  lubrication  equation  one  step  before 


applying  the  differential  discretization  is  given  in 
Eq . { 1 A ) . 
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Discretization  of  Eq.  (lA)  requires  finite  difference 
approximations  of  the  first  and  second  derivatives  of 
Q and  H.  At  this  point,  one  can  already  note  that  the 
differential  discretization  will  have  some  deficiencies 
in  handling  cases  where  significant  discontinuities  of 
clearance  and  clearance-slope  appear. 

The  integral  approach  approximates  Reynolds  equation 
in  the  following  way 


f 1_ 

I 2 


6k,  , (JO 
^*^3  Y 

^..,9  H, 
2Qg  Y> 

. . I 
< .i-f 

+ 

;i>  - 

24|) 

-V/i' 

1 

< , 1 

+ 

t '41  - 

- 

v^Q  - 

■ / 

< -Jf  i 

+ 

;y'  <41  - 

_ ^ 

J 

. 1 . 

< * ? . f 

= 0 

(2A) 

where  A 
A 


y 

X 


6pUyf/(p^h'^.^) 


with  this  method  of  discretization,  only  first  order 
derivatives  of  Q and  H are  needed,  and  here,  at  the  mid- 
points of  the  grid  lines.  It  seems  that  the  integral 
discretization  approximates  the  Reynolds  equation  in  a 
way  such  that  it  provides  the  needed  conservation  of  mass. 

After  the  linearization  of  equations  (lA)  and  (2A) 
each  method  provides  a system  of  linear  equations  that  is 
solved  by  the  so-called  Column  Method.  Both  the  lineari- 
zation and  the  solution  of  the  linear  system  is  described 
in  Ref . [ 9 ] . 

The  rules  that  were  followed  during  the  distribution 
of  grid  points  and  the  variable  grid  spacings  are  outlined 
in  Fig.  3A.  To  achieve  convergence  and  to  eliminate 
numerical  instabilities  a sufficient  number  of  grid  points 
were  needed  to  handle  sharp  pressure  changes  in  boundary 
layers.  In  addition,  since  discontinuities  of  clearance 
slope  exist  in  the  bearing  tested,  it  was  found  to  be 
important  to  position  grid  lines  at  those  points. 

The  two  different  numerical  approximations,  the 
differential  and  the  integral  discretization,  were  first 
applied  to  a tapered  slider  bearing  and  the  pressure 
profile  resulting  from  each  method  agreed  with  each 
other  to  the  third  significant  figure.  Further  testing 
of  these  numerical  methods,  applying  them  to  the  bearing 
under  test,  showed  the  expected  superiority  of  the 
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integral  over  the  differential  discretization.  Results 
for  this  comparison,  given  in  Table  lA,  indicate  a much 
stronger  grid-size  dependence  on  convergence  for  the 
differential  approach. 

•At  this  point,  a final  testing  of  the  numerical 
procedures  was  required.  The  exact,  closed-form  solution 
for  the  infinite-wide  tapered,  tapered-f lat , and  tapered- 
tapered  slider  bearing  with  slip  boundary  condition  was 
obtained.  This  analytical  solution,  outlined  in  Appendix 
B,  was  found  to  be  essential  in  order  to  verify  the 
accuracy  of  the  numerical  solution.  Having  this  exact 
solution  facilitated  a comparison  between  the  two 
numerical  methods  and  the  exact  solution.  Results  giving 
values  of  maximum  pressure  obtained  are  given  in  Tables 
2A,  and  3A.  Fig.  4A  describes  the  pressure  profile 
obtained  using  the  exact  solution  against  results  obtain- 
ed using  the  numerical  procedure  for  a bearing  with 
similar  dimensions  to  those  used  in  the  experimental  work. 
Fig.  4A  also  verifies  the  numerical  solution  for  handling 
both  boundary  layers  and  slip  boundary  conditions. 


_ . 


APPENDIX  B 


Ex^t  Solution  of  the  1-D  Reynolds  Equation 
with  Slip  Boundary  Condition  for  Tapered, 

T a t JC rod-Flat,  and  Tapered-Tapored  Slider  Bearings 


The  solution  for  the  modified  Reynolds  equation  tor 
a tapered,  infinitely  wide  slider  bearing  was  partial  1 
given  by  Burgdorfer  [10] . Because  of  some  errors  in  tin 
solution  in  Ref.  [10],  the  solution  is  repeated  here. 
Slider  bearings,  for  which  the  solution  is  outlined,  can 
have  any  combination  of  tapered  and/or  flat  surfaces. 
Separate  solutions  for  tapered  and  flat  sliders  are  givin 
and  the  conditions  of  pressure  and  mass  flow  for  a com- 
bined bearing  are  shown.  With  the  solutions  described 
here,  one  can  obtain  rhe  bearing  performance  for  any 
bearing  that  has  a combination  of  tapered  and  flat 
surf  aces . 

The  non-dimensional  Reynolds  equation  with  slip 
boundary  condition  for  a slider  bearing  is  given  as 


d 

dX 


'L’JI 

A dX^^ 


Pll 


0 


111-' 


whore  Y 


x/f , H = h/h, , k = X /h 
^ SI 
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A.  Solution  for  a Tapered  Bearing  - Fig.  IB 


Integrating  equation  (IB)  gives 


- \PH  = -K^A 


(2B) 


changing  variables  from  X to  H and  integrating  again  gives 
an  equation  in  the  form 


0^  - (A*  - 6k)9  + A*Kj^  = 


(3B) 


in  which  Kj  and  C are  the  two  constants  of  integration 
and  4' (9)  assumes  the  following  values: 


(0)  = EXP 


★ 1 

A + 6k  ,0 


- ('  /2  -3k) 
/K/2 


(4B) 


for  K > 0 


f2(0)  = EXP 


A + 3k 
l0  - A*/2  + 3k 


(SB) 


for  K = 0 


'('3(0)  = 


0 - _A*/1_L_3 1^+  *^-K/2 
0 - \ */2~+  3k  - /^/2 


A t 6k 


(6B) 


for  K < 0 


£1 
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where  K = 4A *K j - (A *~6k) ^ 

A*=  A/(H,-1) 

0 = PH 

and  Hq=  hfl/h) 

B.  Solution  for  a flat  slider  bearing  - Fig.  2B 

Equation  (IB)  for  a flat  slider  bearing  takes  the 
following  form 


dx'  P Mx  ^ ^ 


(7B] 


Integration  of  Eq . (7B)  gives 


P(1  + 6k/P)^  - AP  = -C^/\ 


(8B) 


Integration  of  Eq . (8B)  gives 


+ (C^/A^  + 6k)Log{P  - C^/A^)  = hX  + 


(9B) 


where  c and  c,  are  the  integration  constants. 
1 2 


C.  Solution  for  a Tapored-Flat  Slider  Bearing  - 
Fig.  3B 


This  solution  for  the  tapered-flat  slider  bearing 
will  demonstrate  how  to  combine  the  solutions  for  the 
tapered  bearing  and  the  solution  for  a flat  bearing.  Thu 
bearing  in  Fig.  3B  is  divided  into  two  regions  indicated 
by  subscript  I and  J.  Solutions  for  the  pressure  in 
those  regions  were  given  in  the  two  previous  sections. 

Eqs . (3B)  and  (9B)  together  indicate  the  pressure.  These 

equations  carry  constants  C^,  C^,  Kj,  and  , to  be 
satisfied  by  boundary  conditions  at  the  leading  and  trail- 
ing edge  of  the  bearing  and  by  conditions  on  pressure 
continuity  and  conservation  of  mass  between  sections  one 
and  two  of  the  bearing. 

The  condition  on  mass  flow  gives 


r 


dP 

[^i"’  ^ p7iT> 


I*'  dx 


I 


6nu  „ 

^2Pa  ^ 


fp  (1  + ^ p 

1 


= 0 


= 0 


(lOB) 


Using  the  condition  on  pressure 


Pi  (Xi=fr)  = 


(1  IB) 


e-  j 
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reduces  Eq.  (lOB)  to  the  following 


■ "/''I 


(12B) 


Finally,  the  four  conditions  needed  to  solve  for  the 
four  constants  in  Eqs.  (3B)  and  (9B)  are  given 


Pf  Pa 

" Pa 

Pl(xi=ij)  = Pj(xj=0) 

Kf  = (13B) 

The  set  of  algebraic  equations  obtained  using  Eqs. 
(3B),  (9B)  and  (13B)  are  transcendental.  Careful  observ- 

ation and  study  of  these  equations  avoids  numerical 
round-off  type  errors  that  can  alter  solutions  signifi- 
cantly, Listing  of  a computer  program  that  solves  Eqs. 
(3B)  and  (6B)  are  given  in  Appendix  C. 


APPENDIX  C 


L i s t i n ‘j  oF  Computer  Pro(;ram  Solving  the  Exact  1-D  Reyno  i • i 
Ec^uation  _wi_th  Slip  Boundary  Condition 


C 

C 

C 

C 

C 

C 

C 

c 

r. 

r 

r. 

c 

c 

c 

f, 

r 

L 

r, 

r. 


TAPtPff)  <^.UDER  HEAR  I NT  ( K . L . 0 . 0 ) 
nij]  PUT  PROV  IPFS  K1  ,K2,  P 

IF  PUT  FLPW  CHART 

RFAD/nATA/Hl,  XI  ,N1  ,P 

<r<f  4<i«ri^4c  *4(4c*  1^4*4  4**#<ri:E*yt 

HI  =U  f API  Ml  E nur  Cl  F AP  AN(  F 1/  ( TR  AL  ING  FCCF  CUARANU) 

R1  =1  Mf)l  fCUL  AR  REAM  FREE  PATH  AT  ARF^.  P ) / I CL  F AR  AN  C t \ I 
T R A 1 1 INC,  F nUF  I 

P=(PPFSSUPF  AT  TRAIlP'G  EDGE  )/(  ARFIEN  PRESSLIFE) 

XI  =6.  ONU^U*!  / ( PA-i<HTR’*HrR  ) 

PA=  A-ThlFNT  PFFSSURE  (PS  I A) 

NU=Eltlin  Vl'-ffSMY  ( Lh^'^  rc/l  IN*  I N I ) 

VI  I UC  I r Y ( 1 N/  SI  C» 

I = Bf  AP  IMF  I r NC.TE' 

HTR=  CLEARANCE  A I TRAILING  EDGE  IIM 

*«  f *«♦*  tt  ♦*♦**♦*  ^***<:***ftt^**ft#«4<;»  ^ 

I PPL  IL  I T R I Al  *('  ( A-H,  n-  /) 

REA|*8  ^1 
PFAI*H  X(2) 

PEAI'-A  ZR  F Al  1 , /Rt  Al  ?,  UF  PT  ST 
EXTFRNAI  F 

CHRMPN  M I , ( 1 . X I . R 1 , P,  Zf  T A ,w  ,r  I 1 ,C?  ? 

N AREl  I S7  /D  A I A/  f 1 , XI  , 'A  1 , P , X , C 1 ,C1  1 , C?? 

I P2-  1 

STF  ps^  zu.npo 
Nsrrps=?i 

1 RE  A|,(  S , DA  I A , I N ( I 

c i = xi  / I HI- 1 .no  ) 

F PS2=0. 001 fO 
E TA  = 0. 01  no 
NS  I G- I 2 
N=  1 

nMflx=l  Ml  *MUG.  rC^HMRI  II  /('  I 

EP=(  1/  ''t  . no  *0.0  C*RI  1*M1  /r  !♦  H 1 *H  1 /C  I ♦ F .nO*R  1*H  1 / ( 1 - A . 

1 Rl-Hl 

20  X<  1 ) =(  OMAX  «F  P)  /2.D0 
x(Pt-(  np»Ax*f  p)  /2.P0 
I TRAX=  200 
W R I I ( I G , 0 A 1 A I 


CALL  ZRE  AL  1(  F,  E PS.fPSZ,  FT  A,NS  ir,,N,  >,  I TKAX.IEP  , ID) 

IF(  lO.EO.l  ) GO  TO  2 
XKl=HI-DWAX*XI  1 ) 

WRITFI/S,  no  OPAX.XKl  ,P 
OEM=XK 1-P 

I F (DABS( DFM) .L F.. lD-3 ) GO  TO  10 
IF(  I.EQ.2I  GO  TO  30 
IF(ID2.F0.l)  GO  TO  10 
P2  = P 

P=P*1. OIDO 
XK2  = XK  I 
1 = 2 

GO  TO  20 
30  PTFMP=P 

P=(P*( XK2-XKI) ♦XKl*(P-P2) )/((P-P2)*(XK?-XK1 ) ) 

P2=PTEMP 
XK2=XK  I 
GO  TO  20 

10  WR1TE(6,110)  DPAX,XK1,P 

no  FORMAT!'  • , • OMAX  = ' rOl  5.  6,  *Kl=  • ,01  5 .6,  20X,'  P =', 

1 015.6,/) 

D3  = H1 

n=!Hl-P)/STEPS 
EP=X(  1 ) 

DO  90  I=1,NSTFFS 

H=(C1*C  (Hl-03)  ♦EP  )*6.  DO  1*(D3-Hl  ) ♦(D3-H1  )<'(  D3«H1  ) )/ 

1 (C22*F IH(D3 ) ) 

H=DSQRT( H) 

XX=IH-Hl  )/  ( l.OCO-Hl  ) 

PRFSS=D3/H 

WPITFI6,120)  PRFSSfH,  XX  ,D3 
120  FORMAT ( • • ,701 7.7) 

D3  = D3-n 
50  CONTINUE 
2 WRITE(6,DATA) 

GO  TO  I 
9959  STOP 
END 

FUNCTION  F<EPS,ID) 

IMPLICIT  REAL*e  IA-H,0-Z) 

REAL08  Ml 

COMMON  H1,C1,XL,M1,P,2F  TA,Q,Cll,C22 

ZFTA=nt.D0*Hl/Cni2.O0»*Ml  /C  1- ( 36.  C C*Ml*M  1 . DO*Hl  *Hn 

1 24.nO*Ml*Hl  )/ICl*Cl)  )*4.D0*EPS/C1 
A=ZETA-l .DC 
10=0 

m A.LE.0.000)  GO  TO  10 
A=OABS( A) 

Q=  CSCRT(A) 

H*H1 


ooooooooono 
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F 

F I2>FIHI PI 

Cll*Cl*EPS/(Hl*Hl*FIll 

C22*(IH1-PI*IC l-P-Hl-6. 00*MII ♦C1*EFSI/FI2 
F*Cll-C22 

C WRITE(6«100)  F,EPS,ZETA,Q,C11,C22»FI1,FI2 

RETURN 

10  WRITEI6.100)  EPS, ZETA.EPS 
100  FORMATC  • ,8015.71 
10*1 
RETURN 
END 

FUNCTION  F IH(THETA) 

IMPLICIT  REAL*8  (A-H,0-Z) 

REAL*8  Ml 

COMMON  H1,C1,XI,M1,P,  ZETA,Q,C11,C22 
A»2.D0*THETA/0/Cl 
8—Cl/O/Cl 
C»6.00*M1/C/C1 
A«0ATAN(  ( A4BKC  I 

F1H»DEXP|-IC146.00*MI 1*2. 000/C1/0*A) 

RETURN 
END 

FUNCTION  - ZREALl  FINOS  REAL  ZEROS  CF  A REAL  FUNCTION 
lER  - ERROR  PARAMETER  (OUTPUT  I 
MARMNG  ERROR  > 32  ♦ N 

N = I INDICATES  THAT  CONVERGENCE  WAS  NUT 
OBTAINED  FOR  AT  LEAST  CNF  INITIAL  GUESS 
(MIThIN  ITMAX  ITERATIONS!. 

X(I ) IS  SET  TO  111111.  FCR  THE  t'S  WHERE 
CONVERGENCE  IS  NOT  CBTAINEO.  NOTE  THAT 
THE  ROUTINE  IS  DESIGNED  SO  THAT  A MULTIPLE 
ROOT  WILL  NOT  APPEAR  IN  THE  OUTPUT 
VECTCR  X. 

SUBROUTINE  ZRFALl(F,EPS,EPS2,tTA,NSIG,N,X,  ITMAX.IFR,  If)) 
DIMENSION  X(  1) 

REAL*8  P,P 1,P2 ,X0,X1, X2 ,RT,FRT ,FPRT,n,CD,DI ,H,BI ,OEN, 

• ON,OM,TEM,X, F ,EPS,FPS2, ETA,CIGT  ,TFN,CNE,ZERC, 

• P9,P11,HALF,PP1,F4 

DATA  TEM,0NE,ZER0,P9, PI  1, HALF  , PP 1 , F A/ 10 . CO, I . DO , 

• 0.00,.9D0,1.  100,.50  0, .100,  A. 00/ 

C DATA  TEN, CNE, ZERO, P9,Pll, half, PPl, FA/10. 0,1. 0,0.0, 

lER  » 0 

OIGT  « TEN**(-NSIGI 

P « -ONE 

PI  » ONE 

P2  « ZERO 

H » ZERO 

DO  95  L * 1,N 

JM  « 0 


/ 


ZERO  GO  TO  5 


IF  (XfLI  .EO. 

P * P9*XIL) 

PI  * P11*X(LI 
P2  = XIL) 

5 RT  * P 
GO  TO  65 
10  IF  UK  .NE.  1)  GO  TO  15 
RT  = PI 
XO  = FPRT 
GO  TO  65 

15  IF  «JK  .NE.  2)  GO  TO  20 
RT  = P2 
XI  = FPRT 
GO  TO  65 

20  IF  UK  ,NE.  3)  GU  TC  55 
X2  = FPRT 
0 = -HALF 

IF  (XIL  » .EG.  ZFRO » GO  TO  25 
H =-PPl*X(L» 

GO  TO  30 
25  H = -ONE 

TO  00  = ONF+0 

81  = XOT'OT'^Z-Xl^DO^^a^XZ^IOD^Dl 

OEN  = 81**2  -FA*X2*0*DD*(X0*0-( X1*0C)*X2) 

IF  ( DEN  .LE  . ZEROI  GO  TO  35 
OEN  = DSQRT (CENl 

GO  TO  40 

35  OEN  = ZERO 

AO  ON  = 81  «■  OEN 

OH  = Bl  - DEN 

IF  (nABS(DN»  .LE.  DARSIDHII  GO  TC  A5 
OEN  = ON 

GO  TO  50 
A5  OEN  = OH 

50  IF  (DEN  .EQ.  ZEROI  DEN  = ONE 

CI=-0C*CX2*X21/0FN 
H = 01  * H 
«T  = RT  ♦ H 

IF  lOARSIHI  .LT.  D AB  SI  RT  I *0  IGT  I GO  TO  90 
GO  TO  65 

55  IF  lOABSIFPRTI  .GE.  OA  BS I X2*10  . CO  11  GC  TO  60 
X0=XI 
X1»X2 
X2»FPRT 

n=oi 

GO  TO  30 

60  Cl  = ni  * HALF 
H = H * HALF 
RT  = RT  - H 
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IF  (JK  .LT.  ITHAX)  GO  TO  75 
IER  = 33 

XIL  ) = 11  11  11  .00 
GO  TO  95 
75  FBT=F1RT,IC) 

FPRT  = FRT 

IF  U .LT.  2)  GO  TO  B1 
on  80  I = ?,L 

TFP  = RT  - X(  I-l) 

IF  lOABSITEMI  ,LT  . EPS2 I GO  If  85 
FPRT  = FPRT/TEM 
PC  CONTINUE 

61  CONTINUE 

IF!  I OARSI  FR  T ) .LT  .EPS  ).  AND.  ( OARS  (FPRT  ) ,LT  .EPS)) 
1 GO  TO  90 


GO  TO  10 

05 

RT  = RT  ♦ 

E TA 

JK  = JK  - 

1 

GO  TO  65 

50 

XIL)  = RT 

55 

CONTINUE 

ITMAX  = JK 

IF(  lER.EO.O) 

GC  TU  9005 

90C0 

CONTINUE 

CALL  UERTST 

( I ER ,6HZR  FA  LI  ) 

90C5 

RETURN 

END 

SUBROUTINE  UER T ST ( I ER ,N AM E ) 

INTEGER*2  NAM E ( 3 ) 

INTEGER  WARN, V*ARF, TERM, PRIMR 
DIMENSION  ITYP  (5,<VI  ,I  BI  T(  9) 

EQUIVALENCE  ( I P I T ( 1 ) , WA RN  ) , ( I P I T ( 2 I , W A R F | , ( 1 P I T ( 3 » , T E RM  I 
DATA  ITYP  /'WAPNSMNG  ', 

♦ 'WARN* ,*  ING( • , • Wl  TH  • , • FIX*,')  ', 

* *TFRM* ,*  INAL* ,'  *,  * * , * * , 

♦ 'NON-* ,* CEF I* ,*NFL  *,*  *,*  '/,  j 

* IBIT  / 32,65,123,0/ 

DATA  PFINTR  / 6/  '1 

IER2=1FR 

IF  nfR2  .GE.  WARN)  GC  TO  5 ' 

IER1=5  ' 

GC  TO  20 

5 IF  (IER2  .LT.  TFRM)  GC  Tl  10 
IER1=3 
GO  TO  20 

1C  IF  HER?  .IT.  WARE)  GO  TO  15 
lER  1*2 
GO  TO  20 
15  IFR1*1 

2C  IFR?»I EP2- IHIT ( IFPl ) 

J 

I 


t 


90 


W«ITP  (PR1NTR,?5J  ( ITYPU  , IFRl),I  = l,5»  ,NAMF,  IEP2,  lER 
25  FOPMATC  *♦*  I M S KUERTSTI  ' ,5A  A , , 3A2 , AX  , 12  , 

* • ( IFR  = • , 13, • ) • ) 

RETURN 

ENT 


APPENDIX  D 


Perturbation  Solution  for  the  1-D  Reynolds  Equation  with 
Slip  Boundary  Condition 

1 D_.  Introduction 

The  modified  Reynolds  equation  describing  the  beha- 
vior of  gas  bearings  in  the  "slip  flow  regime"  was  first 
introduced  by  Burgdorfer  [10],  For  bearings  where  the 
clearance  is  on  the  order  of  a few  mean  free  paths  this 
modified  equation  predicts  a reduction  in  load  carrying 
capacity.  Slip  can  occur  because  of  either  extremely 
thin  bearing  films  or  low  ambient  pressure  conditions. 

The  use  of  singular  perturbation  techniques  for  bearings 
operating  with  large  bearing  numbers  and  including  the 
molecular  mean  free  path  effects  is  here  developed. 
Numerous  previous  studies,  by  Gross  and  Zachmanoglou  (24], 
DiPrima  [25],  and  Schmitt  and  DiPrima  [26]  are  concerned 
with  behavior  of  Reynolds  equation  with  large  bearing 
numbers.  Elrod  and  McCabe  [27]  solved  for  the  asymptotic 
behavior  by  linearization  outside  the  boundary  layer  and 
by  using  the  exact  Harrison  solution  for  a linear  film 
thickness  inside  the  boundary  layer.  None  of  these 
studies  attempt  to  investigate  the  asymptotic  behavior 
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of  bearings  with  large  bearing  numbers  while  including 
the  molecular  mean  free  path  terms. 

The  matched  asymptotic  expansion  technique  laid  out 
by  DiPrima  [25]  is  followed  in  this  analysis.  The  solu- 
tion obtained  in  Sec.  2D,  with  the  Knudsen  number  k set 
to  zero,  is  similar  to  DiPrima' s development;  the  outer 
expansion  and  the  combined  inner  and  outer  solution  are 
correct  in  the  outer  region  to  terms  of  0(1/A^)  and  in 
the  boundary  layer  to  terms  of  0(1/A). 

The  solutions  obtained  in  Sec.  2D,  are  applied  first 
to  a tapered  slider  and  then  to  a parabolic  one.  The 
solutions  for  the  tapered  slider  are  compared  with  the 
exact  analytical  solution  obtained  by  Burgdorfer  [10]. 

The  solution  for  the  parabolic  slider  is  compared  with  a 
numerical  solution.  In  the  latter  case,  it  was  worth- 
while to  use  the  asymptotic  solution  as  the  starting 
solution  for  the  iterative  numerical  scheme. 

2D.  The  Asymptotic  Expansion 

The  use  of  the  modified  Reynolds  equation  to  account 
for  slip  at  the  boundary  for  a slider  bearing,  as 
described  in  Fig.  ID,  was  first  introduced  by  Burgdorfer. 
In  one  dimension,  Burgdorfer  obtained  the  following  non- 
dimensional  nonlinear  ordinary  differential  equation 
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d 

dX 


H ^ P ( 1 


+ 


PH'  dX 


dPH 
'd  X 


(ID) 


where  the  pressure  made  nondimensional  with  respect  to 

the  ambient  pressure  p ,H  is  the  continuously  differen- 

tiable  height  nondimensionalized  with  respect  to  the 

minimum  clearance  h and  X is  the  spatial  variable  in  the 

m 

direction  of  motion  of  the  running  surface  nondimension- 
alized with  respect  to  bearing  length  t.  In  equation 
(ID)  A = 6uU£/(p  hM  is  the  bearing  number  and  k = X^/h, 
is  the  Knudsen  number. 


X=0  X=l. 

Fig.  ID 


V 


The  following  boundary  conditions  for  Eq.  (ID) 
satisfy  the  physical  requirement  of  continuity  of  pressure 


P{X=0)  = 1.0 
P(x=l)  = 1.0 


(2D) 


For  large  speeds  and  small  minimum  clearances,  i.e.  for 
large  compressibility  numbers,  /V , P is  assumed  to  have  an 
expansion  of  the  form 


P°(X;k,A)  = P' 


(3D) 


where  the  superscript  ° is  used  to  indicate  the  pressure 
solution  out  of  the  boundary  layer  with  thickness  of 
order  1/A.  Substitute  (3D)  into  (ID)  and  get 


(p“o4pi  + - • • + - 


.)+H^6k^(p“+ip«  + ...)] 


= ^|x(HP«+H^P°+H;^P»  + ...) 


(4D) 


Eg.  (4D)  describes  the  pressure  behavior  in  the 
region  outside  the  boundary  layer  at  X = 1.0  for  large 
bearing  numbers.  Hence,  the  appropriate  bounday 
condition  for  P®  is 


P“ (X=0.0)  =1.0 


(5D) 


Now,  Eq.  (4D)  is  separated  into  equations  having  terms  of 
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equal  order  with  respect  to  l/^  and  it  becomes  a set  of 
equations  of  the  form, 


0 = 


dHP„ 

air 


j dp“  dP° 

dX^^^^od"3T  “^^’^dir 


) = 


dHP, 
d X 


dP 


dP 


dP^ 


d{HP“) 


dX^H^Po^X-  ^ ^ 


(6Da) 

(6Db) 

(6Dc) 


Here  the  boundary  conditions  (5D)  are 
satisfied  by  the  zeroth  order  term  of  the  pressure,  i.e. 
Pq  (X  = 0)  = 1.0,  and  for  the  higher  order  terms  we  apply 
P°  (X  = 0)  = 0.0  with  n being  larger  than  or  equal  to  1. 
With  those  boundary  conditions  in  mind,  solving  equations 
(6Da) , (6Db) , and  (6Dc) , we  get 


0 

P"  = — 

0 H 

H?  6kH„ 

P?  = lf(H;  - H’)  + - HI 


H, 


(7Da) 

(7Db) 


P”  = qf(H  + 6k) [ (H  + 6k) (H*  + H"H  - H"H  - H'H' 

2 H 0 0 0 0 0 


- H„H'  (H^  - H'  ) ] (7Dc) 

ho 

where  Hj  = — is  the  nondimensional  bearing  clearance  at 
the  leading  edge. 

These  equations,  representing  the  solution  for  the 
pressure  away  from  the  boundary  layer  show  the  effect  of 
the  slip  boundary  condition  to  be  on  the  order  of  1/A , 
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For  the  solution  of  Eq.  (ID) , valid  in  the  trailing 
edge  boundary  layer,  the  stretched  coordinate  ^ = A (1-X) 
is  introduced.  This  change  of  variable  follows  the  argu- 
ment by  DiPrima  [25] . Eq.  (ID)  with  the  stretched 
coordinate  ^ takes  the  form 


d^ 


[ (H^P^  + H^6k)^  1 


d (P^H) 
d ^ 


(8D) 


where  the  superscript  i indicates  the  pressure  inside  the 
boundary  layer. 

Taylor  series  of  H(l-^)  about  C = 0 (X  = 1)  gives 

A 


H(l-|.)  = H,  - ih’  + ^(5/A)'HV  +... 

and  the  pressure  expansion  gives 


(9D) 


= P"  ^ M * ■■■ 


(lOD) 


Substituting  expansions  for  H and  into  eq.  (8D) 
results  in 


= -dT 


•dP 


dp;^ 


dPoH, 


(llDa) 


. -dPj  dPt  .dPo 
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dPf 


dP| 


6k(HtgT-  ■ 2CH,H;6kj-^)  ] 


d 


[P^Hl 


PfH,] 


(llDb) 


where  Eqs.  (llDa)  and  (llDb)  are  the  zeroth  and  the  first 
order  terms  with  respect  to  1/A.  According  to  Eq.  (2D) 
one  of  the  two  needed  boundary  conditions  for  Eqs.  (llDa) 
and  (llDb)  take  the  following  form 


P^(C=0)  = 1.0 

(C=0)  = 0.0  n = l,2,... 


(12D) 


The  secondary  boundary  condition  needed  for  those  second 
order  differential  equations  is  obtained  by  matching  this 
boundary  layer  solution  with  the  outer  solution. 

The  solution  of  the  second  order  differential 
equation  (llDa)  for  the  zeroth  order  inner  expansion 
satisfying  the  boundary  condition  (12D),  gives 


- 1 + 


c - P^ 


(c  + 6k)Log^  _ j 


= -K 


(13D) 


where  c is  the  second  condition  to  be  obtained  by  ordinary 
matching  between  the  outer  solution,  P°  and  the  inner 
solution,  P^ . Following  the  matching  principle  outlined 
in  Chaprters  IV  and  V in  [28]  the  limit  of  Pj  as  x 1 
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is  set  equal  to  the  limit  of  Pq  as  ^ Hence,  the 

missing  condition  is  supplied  by  matching  the  two 
expansions.  This  matching  gives 


(14D) 


c = H 


To  obtain  a composite  expansion,  the  method  of 
additive  composition  is  used.  The  sum  of  the  inner 
and  outer  expansions  is  corrected  by  subtracting  the  part 
that  the  two  expansions  have  in  common.  The  composite 
representation  of  the  two  terms  of  the  outer  expansion  is 
given  by 


P = 


+[p;(a  - ^ (H„  + 6k)  (h;  - H')  (15D) 


where  the  expression  for  Pq is  given  by  Eqs.  (13D) 

★ 

and  (14D).  This  composite  expansion  has  at  least  the 
accuracy  of  each  of  its  constituents. 

Thus,  expression  (15D)  is  in  error  by  no  more  than 
0(1/A^)  away  from  the  boundary  layer  and  0(1/A)  near  the 
boundary  layer.  Although  the  pressure  P^  appears  impli- 
citly in  (13Db)f  the  pressure  evaluation  using  (15D)  is 


see  (25)  and  [28]  for  a more  complete  discussion  on 
expansions  and  matching. 
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straightforward.  After  observing  that  approaches 
exponential  behavior  with  increasing  this  behavior 
can  be  applied  to  obtain  X explicitly  for  any  P<.  Hence 
Eg.  (15D)  can  be  solved  explicitly  by  choosing  a P^ 
between  1 and  Hq,  obtaining  X from  (13D)  using  these  values 
in  (15D) . Using  Eq . (15D)  we  obtain  the  following  form 

for  the  nondimensional  load 


W (a,e  ,k  ) 
H 


/g(P-l)dX  = 


;; (-^  - i)dx  -f,  (H„  - p; (c) )dx  + 


H (H  +6k)  H'-H*  (16D) 

“ ‘ -L- dX 


H 


It  is  useful  to  subdivide  the  load  W in  a way 
similar  to  that  of  DiPrima  [25];  this  subdivision  allows 
direct  comparison  between  the  expansion  for  k equal  to 
zero  [25]  and  k different  from  zero.  Hence, 

W'  = w'  - iw'  + iw  (17D) 

” A ® A 

where  the  first  term  in  Eq.  (16D)  , is  the  limiting 

nondimensional  load  as  A W',  the  last  term  in 

Eg.  (16D) , is  the  value  of  the  load  correction  outside 
of  the  boundary  layer  due  to  terms  on  the  order  of  1/A 
and  W^,  the  second  integral  in  (16D) , is  the  correction 
to  due  to  the  boundary  layer. 
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Using  the  fact  that  this  solution  is  valid  for  a 
large  bearing  number,  it  can  be  shown  for  a parabolic 
slider  that* 


w;  = [Ho  - pf(C)  ]dX  = 


e + + 6ke  + 0 (e  ) (18D) 


where  the  nondimensional  bearing  clearance  H has  the  form 


H = 1 + e [aX^  - (1  + a) X +1] 


(19D) 


Eqs.  (13D),  (15D)  and  (16D)  are  easily  computed  by 

programmable  hand  calculator  or  by  a simple  computer 
program.  Two  simple  examples  are  given  here  in  Sec.  3D. 


Follow  Appendix  B of  [25l  with  some  modification 
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3D.  Examples 

a.  Tapered  slider 

For  this  bearing,  the  analytical  solution  given  by 
Burgdorfer  [10]  was  used  to  check  the  expansion  solution. 
For  this  case,  since  Hj  = Hj  = H' , the  solution  of  (15D) 
and  (16D)  is  significantly  simplified.  The  comparison  for 
the  pressure  profile  given  in  Figs.  2D  and  3D  for  different 
bearing  numbers  shows  a good  correlation  with  the 
analytical  solution.  In  Fig.  4D,  also  a solution  for 
tapered  sliders,  it  is  shown  that  the  reduction  in 
pressure  due  to  slip  vanishes  with  increasing  bearing 
numbers . 


b.  Parabolic  slider 


A slider  with  parc±>olic  film  thickness  H of  the  form 
given  in  Eq.  (19D)  was  investigated.  Sample  calculations 
for  various  values  of  a,  k,  and  the  bearing  number  A with 
e = 2 are  given  in  Table  ID.  A more  detailed  description 
of  the  calculations  are  given  by  DiPrima  [26]  for  the 
case  of  k = 0.0.  The  additional  terms  due  to  mean  free 
path  effects  present  in  Eqs.  (15D)  and  (16D)  did  not  make 
the  calculations  any  more  difficult.  The  load  carrying 
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Pressure  Distribution,  Analytical  vs  Asymototic 


sure  Distribution,  Analytical  vs  Asymptotic 


Figure  4D  Load  vs.  Beciring  Number  for  e = 2.0  and  Various 


= 50 


1000000 


6580 


8422 


1.0265 


,6051 


7095 


7617 


.8139 


. 5220 


6081 


5279 


,6281 


5212 


4336 


.4042 


1.0265 


5258 


4760 


.6481 

.5145 

.4042 

Table  ID  Sample  Calculations  of  W for  Various 
Values  of  a and  ^ (k=0.,l/6) 


■ 
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capacity  does  not  vary  while  changing  the  parameter  "a" 
of  the  parabolic  bearing  profile.  This  behavior  can  be 
applied  in  an  optimal  design  for  constant  loading  of  a 
flat  slider  bearing  where  manufacturing  tolerances  cause 
some  crown  to  exist. 

4D.  Conclusion 


A simple  asymptotic  solution  for  a complicated 
nonlinear  differential  equation  has  many  uses.  The  most 
involved  numerical  solution  of  a nonlinear  differential 
equation  can  be  easily  checked  when  there  exists  some 
knowledge  of  its  asymptotic  behavior.  Here,  such  an 
asymptotic  solution  was  presented.  For  an  infinite 
slider  bearing,  this  solution  predicts  the  limiting  load 
carrying  capacity  with  molecular  mean  free  path  effects 
included . 

The  systematic  application  of  the  asymptotic  expan- 
sion arrives  at  expressions  such  as  Eqs.  (13D) , (15D), 

and  (16D),  which  are  easily  solved  with  a hand-held 
programmable  calculator. 

The  only  restriction  to  this  asymptotic  solution, 
exact  within  0(1/A),  is  the  existence  of  the  first 
derivative  of  the  bearing  clearance.  Following  Schmitt 
and  DiPrima's  work  [26)  and  the  development  that  was 


r-  - 


■* 


given  in  this  paper,  the  asymptotic  method  with  mean  free 
path  effects  could  be  expanded  to  be  applicable  to 
bearings  with  a discontinuity  in  the  film  profile. 

It  should  also  be  noted  that  no  check  of  this 
asymptotic  solution  was  m^de  for  values  of  k greater 
than  1/3.  Although,  purely  from  the  mathematical  point 
of  view,  k could  have  any  value,  the  actual  physics  of 
the  flow  dictate  that  the  modified  Reynolds  equation 
(ID)  is  only  valid  for  small  values  of  k. 
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VII.  NOMENCLATURE 


a 

c 


P 


Pa 

Q 

R 


®1'®2 


= crovfu  parameter  for  parabolic  profile  slider 
= matching  constant 
= integration  constants 
= bearing  clearance 

= bearing  clearance  at  leading  edge 
= bearing  clearance  at  trailing  edge 
= ramp  height;  (Fig.  9) 

= crown  height  of  bearing  skate;  (Fig.  9) 

= nondimensional  bearing  clearance; 

(H=h/hj^,  Ho=hQ/h^,  H^=h^/hj^) 

= ambient  Knudsen  number;  (k=X  /h  ) 

a 1 

= 4A*K^- (A*-6k) ^ 

= integration  constants 
= total  length  of  slider 
= tapered  length 

= distance  from  trailing  edge  to  pivot  point; 
(Fig.  9) 

= nondimensional  pressure  with  respect  to  the 
ambient  pressure;  (P=p/p  ) 

a 

= ambient  pressure 

= (PH)2 

= radial  position  of  slider  center 
= width  of  skates;  (Fig.  9) 


/ 

■A 
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U 


w 


w' 

A. 

x,y 


X 


Y 


Ax . 

4. 


Ayy 


G 

0 

A 

A 

A 

A 


U 


o 


Linear  velocity  of  disk  under  slider 

linear  velocity  in  the  x and  y directions 

width  of  slider;  (Fig.  9 ) 

nondimens ional  load,  W = load/(Pgfw) 

unit  vectors  in  the  x and  y directions 

nondimensional  spatial  coordinate  in  the  direction 

of  motion 

nondimensional  spatial  coordinate  in  the  direction 
perpendicular  to  motion 

grid  spacing  in  the  x direction  between  points 
i and  -t-1 

grid  spacing  in  the  y direction  between  points 
j and  j-1 

parameter  for  parabolic  slider 
PH 

wave  length  of  monochromatic  light  source 

molecular  mean  free  path  at  ambient  pressure 

bearing  number  = 6 jjU£/ (p,h,  ^ ) 

a i 

A/(Ho-1) 

2 

bearing  number  = 6uU  £/(p^hi  ) 

X 3 X 

bearing  number  = 6yUy£/ (pghj^^) 
viscosity  of  lubricant 
streched  boundary  layer  coordinate 
surface  accommodation  coefficient 


/I 


i 


transcendential  function  of  0 


Differential  Operator 


V 


Superscripts 

0 = refers  to  the  outer  solution 

1 = refers  to  the  inner  solution 

' = d/dX 


Subscripts 


J 


ambient  conditions 
order  of  perturbation 
x-y  index  of  grid  points 

refers  to  values  for  the  tapered  region  of  the 
tapered-flat  slider  bearing 

refers  to  values  for  the  flat  region  of  the 
tapered-flat  slider  bearing 


Ill 
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The  experimental  investigation  described  here  involves  the  highly  accurate 
measurement  of  bearing  clearances  on  the  order  of  10  microinches  in  self-acting 
pivoted  narrow-slider  gas  bearings.  The  experimental  measurements  are  based  on 
light  interferometry  using  a variable-wavelength  pulsed  dye  laser  and  a CW  HeNe 
laser  as  monochromatic  sources.  The  light  interference  in  the  gas  bearing  is 
obtained  by  flying  the  slider  on  a very  precise  optically  flat  quartz  disk  through 
which  the  light  beam  is  transmitted. 

The  combined  effect  of  high  Knudsen  numbers  and  surface  irregularities  on  the 
flying  height  of  narrow  gas  bearings  is  observed  by  varying  the  load  on  the  bear- 
ing and  the  ambient  molecular  mean  free  path.  The  experimentally  measured  bearing 
clearances  are  compared  quantitatively  with  rather  accurate  theoretical  predictions 
obtained  by  numerical  solution  of  Reynolds  differential  equation  for  compressible 
fluids  with  slip  boundary  conditions. 

The  result  of  this  study  indicates  that,  as  clearances  in  narrow  gas  bearings 
get  progressively  smaller,  while  the  Knudsen  number  increases  to  values  beyond  0.1, 
the  theoretical  model  fails  to  predict  the  bearing  behavior.  It  is  also  argued 
that  this  failure  is  because  of  the  weakness  of  the  continuum  model. 
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